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Abstract 

We consider an infinite one-dimensional anisotropic XY spin chain with a nearest-neighbor time- 
dependent Heisenberg coupling J(t) between the spins in presence of a time-dependent magnetic 
field h(t). We discuss a general solution for the system and present an exact solution for particular 
choice of J and h of practical interest. We investigate the dynamics of entanglement for different 
degrees of anisotropy of the system and at both zero and finite temperatures. We find that the 
time evolution of entanglement in the system shows non-ergodic and critical behavior at zero and 
finite temperatures and different degrees of anisotropy. The asymptotic behavior of entanglement 
at the infinite time limit at zero temperature and constant J and h depends only the parameter 
A = J/h rather than the individual values of J and h for all degrees of anisotropy but changes for 
nonzero temperature. Furthermore, the asymptotic behavior is very sensitive to the initial values 
of J and h and for particular choices we may create finite asymptotic entanglement regardless of 
the final values of J and h. The persistence of quantum effects in the system as it evolves and as 
the temperature is raised is studied by monitoring the entanglement. We find that the quantum 
effects dominate within certain regions of the kT-X space that vary significantly depending on the 
degree of the anisotropy of the system. Particularly, the quantum effects in the Ising model case 
persist in the vicinity of both its critical phase transition point and zero temperature as it evolves 
in time. Moreover, the interplay between the different system parameters to tune and control the 
entanglement evolution is explored. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 
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I. INTRODUCTION 



Quantum entanglement represents one of the corner stones of the quantum mechanics 
theory and is of fundamental interest in modern physics 1|. In the early days of the theory, 
the notion of entanglement was first noted and introduced by Einstein, Podolsky, and Rosen 
as a paradox in the formalism of the quantum theory Q. Nowadays entanglement is treated 
as a well-established concept and experimentally verified phenomenon in modern physics. 
Quantum entanglement is a nonlocal correlation between two (or more) quantum systems 
such that the description of their states has to be done with reference to each other even if 
they are spatially well separated. Understanding and quantifying entanglement may provide 
an answer for many questions regarding the behavior of complex quantum systems. Par- 
ticularly, entanglement is considered as the physical property responsible for the long-range 
quantum correlations accompanying a quantum phase transition in many-body systems at 
zero temperature 0-5]. Particular fields where entanglement plays a crucial role are quan- 
tum teleportation, quantum cryptography, and quantum computing, where it is considered 
as the physical basis for manipulating linear superpositions of quantum states to implement 
the different proposed quantum computing algorithms 

Different physical systems have been proposed as reliable candidates for the underlying 
technology of quantum computing and quantum information processing js-16]. The basic 
idea in each one of these systems is to define certain quantum degree of freedom to serve as 
a qubit, such as the charge, orbital, or spin angular momentum. This is usually followed by 
finding a controllable mechanism to form an entanglement between a two-qubit system in 
such a way to produce a fundamental quantum computing gate such as an exclusive Boolean 
XOR. In addition, we have to be able to coherently manipulate such an entangled state to 
provide an efficient computational process. Such coherent manipulation of entangled states 
has been observed in different systems such as isolated trapped ions [l^] and superconducting 
junctions 18[. The coherent control of a two-electron spin state in a coupled quantum dot 
was achieved experimentally, in which the coupling mechanism is the Heisenberg exchange 



interaction between the electron spins 19|-|2l|. 



The obvious demand for a controllable mechanism led to one of the most interesting 
proposals in that regard which is to introduce a time-dependent exchange interaction be- 
tween the two valence spins on a doubled quantum dot system as the coupling mechanism 
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23|. The coupling can be pulsed over definite intervals resulting a swap gate which 



can be achieved by raising and lowering the potential barrier between the two dots through 
controllable gate voltage. The ground state of the two-coupled electrons is a spin singlet, 
which is a highly entangled spin state. 

The interacting Heisenberg spin chain model represents a very reliable model for con- 
structing quantum computing schemes in different solid state systems and a very rich model 
for studying the novel physics of localized spin systems. This spin chain can be experimen- 
tally realized, for instance, as a one-dimensional chain of coupled nano quantum dots. 

There has been many studies focusing on the entanglement at zero and finite temperature 
for isotropic and anisotropic Heisenberg spin chains in presence and absence of an external 



magnetic field 



24J-|33| . Particularly, the dynamics of thermal entanglement has been studied 



in an XY spin chain considering a constant nearest-neighbor exchange interaction, in the 
presence of a time- varying; magnetic field represented by a step, exponential, and sinusoidal 
functions of time [34J, 135 ]. 

Recently, the dynamics of entanglement in a one-dimensional Ising spin chain at zer o 
temperature was investigated numerically where the number of spins was seven at most 



The generation and transportation of the entanglement through the chain under the effect of 
an external magnetic field and irradiated by a weak resonant field were studied. It was shown 
that the remote entanglement between the spins is generated and transported although only 
nearest-neighbor coupling was considered. Later, the anisotropic XY model for a small 
number of spins, with a time-dependent nearest-neighbor coupling at zero temperature was 
studied, too js^j. The time-dependent spin-spin coupling was represented by a dc part and a 
sinusoidal ac part. It was found that there is an entanglement resonance through the chain 
whenever the ac coupling frequency matches the Zeeman splitting. 

In this work, we investigate the evolution of quantum entanglement in an infinite one- 
dimensional XY spin chain system coupled through nearest-neighbor interaction under the 
effect of a time- varying magnetic field h(t) at zero and finite temperature. We consider a 
time- dependent nearest-neighbor Heisenberg coupling J(t) between the spins on the chain. 
We discuss a general solution for the problem for any time dependence form of the coupling 
and magnetic field and present an exact solution for a particular case of practical interest, 
namely a step function form for both the coupling and the magnetic field. We focus on 
the dynamics of entanglement between any two spins in the chain and its asymptotic be- 
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havior under the interplay of the time-dependent coupling and magnetic field. Moreover, 
we investigate the persistence of quantum effects in the system as it evolves in time and 
as its temperature increases. We show that the time evolution and asymptotic behavior of 
entanglement for static coupling and magnetic field at zero temperature depends only the 
ratio of the coupling to the magnetic field rather than their individual values but not at 
finite temperatures. The entanglement was found to be very sensitive to the initial values 
of the coupling and the magnetic field and in particular cases they may dictate the asymp- 
totic entanglement regardless of the final values of the parameters. The quantum effects 
were shown to dominate within certain regions of the temperature, coupling, and magnetic 
field space which depend significantly on the degree of anisotropy of the coupling which are 
manifested by the asymptotic behavior of entanglement. 

This article is organized as follows. In Sec. II we present our model and discuss a 
general solution for the the XY spin chain for a general form of the coupling and magnetic 
field and focus on a particular case where the system is exactly solvable. In Sec. Ill we 
evaluate the magnetization and the spin-spin correlation functions of the system and use 
them to evaluate the entanglement. In sees. IV, V, and VI we study the entanglement 
dynamics in the completely anisotropic, partially anisotropic, and isotropic cases of the 
system respectively. We conclude in Sec. VII and discuss future directions. 

II. THE TIME DEPENDENT XY MODEL 

In this section, we present an exact solution for the XY model of a one-dimensional 
lattice with N sites in a time-dependent external magnetic field hit). We consider a time- 
dependent coupling J(t) between the nearest-neighbor spins on the chain. The Hamiltonian 
for such a system is given by 

Jit) N Jit) N N 

H = + 7) E ^X+i - - 7) E ~ E Kt)at , (1) 

z i=i z i=\ i=i 

where <jj's are the Pauli matrices and 7 is the anisotropy parameter. For simplicity, we will 
consider Ti = 1 throughout this article. Introducing the raising and lowering operators aj, aj 

<4 = \(of + io1), a^\{^-ia\). (2) 
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Hence, Pauli matrices can be written as follows: 



T t V Z <r\ f T 

<jj = a- + en, of = : — , cjj = 2a\ai — 1 . 



(3) 



Following the standard procedure to treat the Hamiltonian ([I]), we transform the Pauli spin 
operators into fermionic creation and annihilation operators b\, hi 38] 



i-l 



i-1 



i\ = b\ exp(i7r ^2 bjbj), a, = exp(— iir b^bj)bi 

3=1 3=1 



(4) 



then, applying a Fourier transformation, we obtain 



b] 



N/2 



N/2 



_L_ e ij<t}p c\ h = — — V 



(5) 



p=-N/2 



where P = As a result, the Hamiltonian can be written as 

N/2 

H = H p , 
P =i 



(6) 



with H p given by 



Hp = a p (t)[c ] p Cp + cLpC-p] + iJ(t)5 p [clcL p + c p c_ p ] + 2/i(t) 



(7) 



where a p {t) = — 2J(t) cos0 p — 2h(t) and 5 P = 27sin0 p . 

The decomposition of the Hamiltonian was only possible because [Hi,H m ] = 0, for 
I, m — 0, 1, 2, . . . , N/2. Thus the Hamiltonian in the 2 7V -dimensional Hilbert space has been 
split into N/2 noncommuting sub-Hamiltonians, each in a four- dimensional independent sub- 
space. Writing the matrix representation of H p in the basis {|0) , c^c_ p |0) , c p |0) , c_ p |0)} we 
obtain 



H n 



2h(t) -iJ(t)8 p 

iJ(t)5 p -AJ (t) cos <p p - 2h{t) 

-2J(t)cos0 p 



(8) 



\ -2J(t) cos (j) p J 

Initially the system is assumed to be in a thermal equilibrium state and therefore its 
initial density matrix is given by 



P P (0) 



-j8Hp(0) 



(9) 



where = 1/kT, k is Boltzmann constant, and T is the temperature. Using Eq.Q the 
matrix representation of p p (0) reads 



D 2/3(cos </> p +r[/i(0),J(0)]) 



where 



^11 



2T[h(0),J(0)\ 



c v ■ 

S21 ■ 



/ Cfi Cf 2 o o ^ 

C!i C 2 P 2 o o 

o o cf 3 o 

V o o o ( 4 P 4 ; 



{T[h{0), J(0)] + J(0) cos0 p + /i(0)} e" 4 ^!^ )-^ )] 

+ {r[/i(o), j(o)] - j(o) cos0 p - /i(o)}] , 

^ p j {l_ e -4/3r[Mo),J(o)]| 



4T[fc(0),J(0)] 



■^22 



2T[h(0),J(0)] 



4I>(0),J(0)] 

{T[h(Q), J(0)] - J(0) cos0 p - /i(0)} e-^Moum 
+ {V[h(0), J(0)] + J(0) cos P + /i(0)}l , 



fP _ A 
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-2/?r[h(o),j(o)] 



and 



r[h(t),J(t)} = {[J{t) cos p + /i(t)] 2 + 7 2 J 2 (t) sin 2 <j) p y . 



(10) 



(11) 
(12) 

(13) 



(14) 
(15) 

(16) 



Since the Hamiltonian is decomposable we can find the density matrix at any time t, p p (t), 
for the pth subspace by solving the Liouville equation given by 



ip p {t) = [H p {t),p p {t)\ 



(17) 



which gives 



p P (t) = U p (t) Pp (0)Ui(t), (18) 
where U p (t) is the time evolution matrix which can be obtained by solving the equation 



iU P {t) = U p (t)H p (t) . 



(19) 



Since H p is block diagonal U p should take the form 



U p (t) 



Ul U? 



'11 w 12 

jP TjP 

'21 u 22 



ul vz, 



V 



f/ 3 p 3 
U p u j 



(20) 



Fortunately, Eq. f lT9|) may have an exact solution for a time-dependent step function form 
for both exchange coupling and the magnetic field which we adopt in this work. Other 
time-dependent function forms will be considered in a future work where other techniques 
can be applied. The coupling and magnetic field are represented respectively by 



J{t) = Jo + ( Ji - J )9(t) 



(21) 



h(t) = h + (hi - h o )0(t) 



(22) 



where 9(t) is the usual mathematical step function. With this set up, the matrix elements 

—i[Ji cos (pp + hi] sin[2tr(/ii, Ji)] 



of Up were evaluated to be 



jjP^ _ g2it Ji cos 4>v | . 



nKh) +cos[2tr(hi, Ji)]} 

— JiS p sin[2tr(hi, J\)] ■> 

J ' 



; 22 



rrp _ 2itJi COS (j> p ( 

Ul2 ~ e t 21^, * 

r rP _ ^MtJicoB^ f sin[2tr(fe x , Jx)] 



"I 



r(/ii, Ji) 

^33 — ^44 — e 



Consequently, the density matrix takes the form 



p p (t) = e 



2/3J o cos0 p +2/3r(/i o ,Jo) 



( f? n p\ 2 ^ 

P21 P22 

p p 33 

V ff M J 



(23) 

(24) 

(25) 

(26) 
(27) 



(28) 
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where 



{MJoh - Ji/io]<jJsin 2 [2tr(/n, Ji)] 



4r(/i , Jo)r 2 (/n, Ji) 

+ 2r 2 (/i 1 , J 1 )[r(/i , Jo) + J o cos0 p + /^je" 4 ^ ^ 

+ M^ho - J /ii]^sin 2 [2ir(/n, J x )] 

+ 2r 2 (/i!, Ji)[T(/io, J ) - J o cos0 p - /io]} , 



(29) 



^(1 _ e-Wo.Jo)) 



4r(/, , Jo)^^, Ji) l r(/l1 ' Jl) Wl - JlM sin[4tr(/l1 ' Jl] 

+i {j T 2 (h u Jx) + 2( Ji/io - J /ii)( Ji cos0 p + /ii) sin 2 [2*r(/n, Ji)] 

P21 = (P12)* , 



(30) 
(31) 



P2 



1 



{Ji[Jiho - Joh}S 2 p sm 2 [2tT(h u Jx)] 



22 ir(/'o.-/o)r 2 (/',.-/,) 

+ 2r 2 (/n, Ji)[r(/i , Jo) - J o cos0 p - h ]}e~^ h ^ 

+ MJoh - Jx/io]^ sin 2 [2*r(/n, Ji)] 

+ 2r 2 (/n, Ji)[r(/i , Jo) + J o cos0 p + /io]} , 



^33 — T44 



-2/8r(/i ,Jo) 



(32) 



III. SPIN CORRELATION FUNCTIONS AND ENTANGLEMENT EVALUA- 
TION 

In this section we evaluate different magnetic functions of the XY model which we utilize 
afterward to evaluate the spin-spin entanglement in the chain. The first function is the 
magnetization in the z direction which is defined as 

1 N 1 1/N 

M=^US-) = ttT,M p , (33) 
JV 3=1 JV p=i 

where M p = c p c p + cL p c_ p — 1. In terms of the density matrix, it is given by 

/nj ,_ Tr[Mp(t)} 1 l ^ Tr[M pPp (t)} 

{ Z) ~ Tr[ P (t)} Tr[ Pp (t)} ' {6A) 
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which yields 



N/2 



tanh[/3r(/i , J )] 



AN^ 1 T^h 1 ,.J 1 )T(h ,J ) 
2J x {J Q h x - Ji/i )5psm 2 [2tr(/n, J x )] + 4r 2 (/i 1 , Ji)(J o cos0 p + V 

The other functions needed are the spin correlation functions defined by 



QX / QX QX \ ny / ny ny \ QZ / QZ QZ \ 



which can be written in terms of the fermionic operators as follows 



(35) 



Sf,m ~ 7 {BiAi + iB i+ i . . . A m _iB m _iA m ) 



D l,rn 



("I) 



l—m 



{AiBi + iAi + i . . . B m ^iA m -iB m ) 



(36) 
(37) 



where 



Using Wick theorem 



Sf, m = J (AiBiA m B m ) 



Ai = bl + bi, Bi = b\-bi. 



(38) 



(39) 



39) , the expressions f l36|) - fl38|) can be evaluated as Pfaffians of the form 



Fii + i Gn + i ■ ■ ■ Gi 

,771— 1 Fl,m 



Pi+ l.z+1 ' ' ' Pi 



i+l,m-l Ql+l,m 



Pm—X,m—1 Qm—l,n 

F i 

u ± m— 1,17 





/ 



c.v 
°Z,7i 



-1) 



l—m 



-Pf 



o Qz,z+i 

Fi + n + x 



Ql,m—X Pl,m 
" l+l,m— 1 ,m 

Pm— l,m— 1 Cm- l,m 
Pm— l,m 





/ 



(40) 



(41) 
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where 



Pi, i Qi, m Pi 

Fl m Glm 



o p 



(42) 



Fi,m = {BiA m ) , P^ m = (A\B m ) , Qi m = (AiA m ) , G i)in = (BiB m ) , (43) 



and 



1 N/2 ( 

Ql,m = T7^2\2 COs[(m — l)(f) p ] 
p=l ^ 

i(Jih — J Q hi)5pSm[(m — 1)4> P ] sin[4£T(/ii, J x )] tanh[/3r(/i , Jo)] 



r(/n, JOrc^cJo) 



(44) 



l ^ /2 r 

G^m = ^7^1 -2cos[(m - l)(j) p ] 

i(J\ho — Johi)5 p sin[(m — /)0 P ] sin[4tr(/ii, Ji)] tanh[/3r(/i , Jo)] 



r(^, jor^o, Jo) 



(45) 



1 ^ tanh[/3r(/i , Jo)] J u A , , 

= ^§P(fc 1> j 1 )r(^Jo)i OOBKTO -^ 1 

J^Jo/i! - Jx/io]^ sin 2 [2tr(/n, Jx)] + 2r 2 (/i 1; Ji)(J cos0 p + ho] 
+5 p sm[(m - l)(j) p ) 
JoT 2 (h u J x ) + 2(J x h - Joh x )(J\ cos0 p + h x ) sin 2 [2tV(h x , Jx)] 



\m J/,m • 



(46) 



The amount of entanglement between two quantum systems, bipartite entanglement, is 



a monotonic function of what is called the concurrence [40] . The concurrence varies from a 
minimum value of zero to a maximum of one coinciding with the entanglement function range 
and behavior. Therefore, the concurrence itself is considered as a measure of entanglement. 
The concurrence C(t) is defined as 



C(p) = max(0, A a - X b - A c - X d ) 



(47) 
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where the Aj's are the positive square root of the eigenvalues, in a descending order, of the 
matrix R defined by 

(48) 



R = \/VppVp: 

and p is the spin-flipped density matrix given by 



p = {a y ®a y )p*{a 



y yy CTy) 



(49) 



Knowing that p is symmetrical and real due to the symmetries of the Hamiltonian and 
particularly the global phase flip symmetry, there will be only six nonzero distinguished 
matrix elements of p which takes the form 4l| 



(50) 



( Pi,i p M ^ 
p 2 ,2 p2,3 
P2,3 p3,3 
\Pia p 4)4 J 

As a result, the roots of the matrix R come out to be A a = • v /pi,iP4,4 + \pi,a\, A5 = v /p 2j 2P3,3 + 



|P2,3|,A C = ^/pTJpaa - \pi,a\ , and \ d = V /P2,2P3,3 - |P2,3 1 

To find the nonzero matrix elements of p, one can utilize the formula of the expectation 
value of an operator in terms of density matrix (^G^ = Tr(p G)/ Tr(p) along with the 
magnetization Eq. (l34"]) and the spin correlation functions Eq. (!36l) - (l38l) which give 



Pi,i = \Mi + \M z m + Sf >m + , 

p 2 ,2 = \m z { - X -M z m - S z hm + - A , 
1 x 1 z 1 

P3,3 = _ ~ ^l,m + 7 ) 



1 



P4,4 



-M? 



1 



Ml 



2 — , 2 - m 
P2,3 = Sf m + Sf 



1 

4 



Pl,4 



Sf 



qy 



(51) 

(52) 

(53) 

(54) 
(55) 
(56) 



12 



IV. TRANSVERSE ISING MODEL 



Considering a completely anisotropic XY model by setting 7 = 1 in the Hamiltonian (JT]), 
we obtain the transverse Ising model Hamiltonian 

N N 

ff=-J(f)£*f*tVi (57) 

-i = l 1 = 1 

Defining a dimensionless coupling parameter A = J/h, the ground state of the Ising model 
is characterized by a quantum phase transition that takes place at A close to the critical 
value A c = 1 [4]. The order parameter is the magnetization (a x ) which differs from zero for 
A > A c and zero otherwise. The ground state of the system is paramagnetic when A — > 
where the spins get aligned in the magnetic field direction, the z direction. For the other 
extreme case when A — > 00 the ground state is ferromagnetic and the spins are all aligned in 
the x direction. The ground state of the Ising model as A — > is a product of individual spin 
states pointing in the z direction, while for A — > 00 is product of spin states pointing in the x 
direction. This means that in both cases the state is minimally entangled. Quantum phase 
transition takes place at zero temperature as the thermal fluctuations destroy the quantum 
correlations in the ground state of the system. The effect of the temperature on entanglement 
near the critical point in the Ising model has been studied in Ref. |4j, where it has been 
shown how the entanglement decays abruptly as the temperature raises; nevertheless, it is 
sustained in the vicinity of the critical point close to kT = 0. In this section we study the 
dynamics of entanglement in the Ising model under the effect of nearest-neighbor coupling 
and external magnetic field where both are considered time-dependent. The number of spins 
N in the system is set to 1000 throughout this study, where testing larger values of N showed 
no effect on the results. 

In Fig. [TJ we explore the dynamics of the nearest-neighbor concurrence C(i,i + 1) at zero 



temperature. In Fig. 1(a) we choose the magnetic field to have a constant value of 1 while 
the coupling parameter takes the value 1 or 0.5 or a step function changing between 0.5 and 



1 (or 1 and 0.5). In Fig. 1(b) we set the coupling parameter to be constant this time with 
a value 1 while the magnetic field can take the values 1 or 0.5 or a step function changing 
between 0.5 and 1 (or 1 and 0.5). As one can see, when the coupling parameter (the magnetic 
field) is a step function, the concurrence reaches a value that is neither its value when J = Jq 
(h = ho) nor J x {h = hi), i.e., concurrence C(i,i + 1) shows a nonergodic behavior. This 
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(a)C 



0.2 



0.1 



(c) ° 



-J =1.J,=1 

J =0.5,^-0.5 
J -1 ,J,.0.5 
- J„»0.5,J.-1 




(b) 0.3 



0.2 



0.1 



___h 


=1,h,-1 


h 


=0.5,^=0.5 


- - h 


=1,h 1 =0.5 


— h 


=0.5,^=1 



, .... .., 



(d) 0.25 



0.20 



0.10 



10 20 

f 




FIG. 1: Dynamics of the nearest-neighbor concurrence C(i,i + 1), where t is in units of Jj -1 , with 7 = 1 
and (a) ho = h\ = 1 for various values of Jo and Ji at fcT = 0; (b) Jq = J\ = 1 for various values of /io 
and /ii at fcT = 0. (c) Dynamics of the magnetization per spin; (d) dynamics of the spin-spin correlation 
function in the z-direction for fixed h = ho = hi = 1 for various values of Jo and J\ at kT = with 7 = 1. 

behavior follows from the nonergodic properties of the magnetization and the spin-spin 
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43). The nonergodic behavior 



correlation functions as reported by previous studies 
of the magnetization and the spin-spin correlation function in the z direction are shown 
in Figs. 



1 c 



and 1(d) The spin-spin correlation functions in the x and y-directions show 
similar behavior. At higher temperatures the nonergodic behavior of the system sustains 
but with reduced magnitude of the asymptotic concurrence (as t — > oo). 

In Fig. [2] we study the behavior of the nearest-neighbor concurrence C(i,i + 1) as a 
function of A for different values of J and h at different temperatures. In this figure we set 



J = Jq = J\ and h — ho — hi. In Fig. 2(a) , we study the zero temperature case, where we 



fix h (or J) and vary J (or h). As one can see, the behavior of C(i,i + 1) depends only on 
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(a) 0.25 
0.2 

^ 0.15 
0.1 
0.05 



(b) 0.25 



- h=0.25 h=1 h=4, kT=1 



5 10 

X 



(C) 0.25 



-J=0.25 J.1 J=4, kT=1 



5 10 

X 



W)o.s 



-h=0.25 h=1 h=4, kT=3 




FIG. 2: C(i,i + 1) as a function of A for h = ho = hi and J = Jq = J\ at (a) kT — with any combination 
of J and h; (b) feT = 1 with h = h x = 0.25, 1,4; (c) fcT = 1 with J = J 1 = 0.25, 1,4; (d) fcT = 3 with 
ho^h^ 0.25,1,4. 

the ratio J/h (i.e., A) rather than on their individual values. As expected C(i,i + 1) starts 
at zero, reaches a maximum value at A ~ A c = 1, and then vanishes for larger values of A. 
Studying entanglement at nonzero temperatures shows that the maximum value of C{i, z+1) 
decreases as the temperature increases. Furthermore, C(i, i + 1) shows a dependence on the 



individual values of J and h, not only their ratio. In Fig. 2(b) one can see that increasing h 



at kT = 1 causes the maximum concurrence to increase and to shift toward smaller values 



of A. In Fig. 2(c), increasing J at kT = 1 causes the maximum concurrence to decrease 



and to move toward larger values of A. Figure 2(d) shows the significant reduction in the 
entanglement as the temperature increases further (kT = 3). 

In Fig. |3] we investigate the dynamics of entanglement as a function of the coupling 
parameter J(t). We plot C(i,i + 1) as a function of time and Ai(= Ji/hi) where we set 
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FIG. 3: (Color online) C(i,i + 1) as a function of Ai and t, in units of J-f , at kT = with 7 = 1, 
h{) = h\ = 1, and (a) Jq = 1; (b) Jo = 5. (c) The asymptotic behavior of C(i,i + 1) as a function of Ai with 
7 = 1, ho = hi = 1, and Jo — 0.5, 1, 2 at kT — 0; (d) dynamics of C(i, i + 1) with 7 = 1 , ho = hi = 1, 
J = 1, 2, and Ji = at fcJ = 0. 



h = h = hi. In Fig. 3(a) we set the parameter values as h — 1, Jq — 1, and kT = 0, i.e., the 



system is initially prepared in a state of maximum entanglement (Ao = A c ). As one can see, 
for zero Ai, the entanglement shows an oscillatory behavior in time, where the spins precess 
about the magnetic field in the z direction, and its magnitude increases as we increase Ai 
until it reaches its maximum value close to A c . As Ai exceeds A c the entanglement decreases 
and eventually vanishes for Ai > 2.5 where in this case J is dominating over h and the spins 
are completely aligned in the x direction. On the other hand, when the system is prepared 
in an initial state with Aq different from A c , the maximum entanglement it can reach is much 



lower than the previous case and appears at lower value of Ai ~ 0.5 as shown in Fig. 3(b) 



where h = 1 and Jq = 5. In Fig. 3(c) we exploit the asymptotic concurrence (as t — > 00) 
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versus Ai for different values of Jo while the magnetic field is set as ho = hi = 1. Clearly, 
the asymptotic value of the concurrence varies significantly depending on the initial values 
of the parameters. For Ao < A c the asymptotic concurrence peaks at Ai ~ 0.5, and for 
A = A c the peak takes place at Ai ~ 1. For higher values of the parameter A > A c two 
peaks show up with a smaller second peak that decreases and shifts toward larger values 
of Ai as Ao increases. For instance, for Ao = 2 the Ai values at the peaks are 0.2 and 1 
with much lower maximum values of concurrence than the previous cases. In Fig. 3(d) we 



examine the dynamical behavior of the concurrence as the coupling parameter is switched 
off, i.e., Jo has a finite value while J\ = 0. As one can see, the oscillation sustains as time 
elapses. Testing different values of h and Jo we find that the oscillation amplitude is largest 
when 0.5 < Ao < 1 while it almost vanishes and loses uniformity when Ao > 2. The period 
of the oscillation decreases as h increases and is independent of Jo- We observe a similar 
behavior for the concurrence when the magnetic field is switched off while setting J = J\. 

In Fig. H] we manifest the dependence of the asymptotic behavior (as t — > oo) of the 
nearest-neighbor concurrence on the magnetic field and coupling parameters ho, hi, Jo, and 



Ji at zero temperature. In Fig. 4(a) we present a three-dimensional plot for the concurrence 
versus Jo and J\ where we set the magnetic field at h = hi = 1. The concurrence starts 
with a zero value for Jo = J\ = and increases abruptly, reaching a maximum value of 
approximately 0.26 at J = J\ ~ 0.88. For J\ > 4, C(i, i + 1) vanishes for all Jo values. It is 
interesting to see that for all initial values Jo > 1 and J\ < 1, the asymptotic concurrence has 
a finite value which decays for higher values of J\. This emphasis that starting with a finite 
coupling Jq and reducing it to a very small value, Ji, leads to persisting entanglement in the 



system. Figure 4(b) shows the asymptotic behavior of the nearest-neighbor concurrence as a 
function of ho and hi, while fixing the coupling parameter at Jo = Ji = 2. The concurrence 
is zero at hi = ho = and increases as both increase until it reaches its maximum value at 
ho = hi ~ 1.8, i.e., close to the A c . If we start with a relatively large magnetic field (say 
ho > 3), C(i,i + 1) will vanish for hi < 1, reach a maximum value at hi ~ 2, and then 
decrease gradually with further increase of hi. However, if we start with a smaller magnetic 
field (1 < ho < 2), C(i, i + 1) will have a maximum value if hi is kept within the same range 
and will vanish when hi is increased. Finally, if we start with a much smaller magnetic 
field (ho << 1), C(i,i + 1) vanishes for hi < 1.5 but increases as hi increases and reaches 
a plateau, i.e., there will be a finite asymptotic concurrence left in the system for a very 
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FIG. 4: (Color online) The asymptotic behavior of C(i, i + 1) as a function of (a) Jo and J\ with 7 = 1 and 
ho — h\ — 1; (b) ho and hi at kT = with 7=1 and Jo = J\ — 2; (c) ho and Jo at fcT = with 7 = 1 and 
hi — Ji = 1; (d) /ii and Ji at kT = with 7 = 1 and ho = Jo = 1, where /iq, far, and Jq are in units of J\. 



small initial magnetic field and large final one. The asymptotic behavior of C(i,i + 1) as 



a function of Jo and ho at kT = is illustrated in Fig. 4(c) We studied this behavior for 
Ji = hi = 0.5, Ji = hi = 1, and Ji = hi = 2 and in all cases we got the same behavior. 
The largest entanglement is reached when J = /io but as Jo differs from h , C(i,i + 1) 
decays in agreement with the physical interpretation discussed previously. We also study 
the asymptotic behavior of C(i,i + 1) as a function of Ji and hi at kT = while setting 
Jo = ho = 0.5, 1, and 2. Again, we get the same behavior in all three cases as shown in 



Fig. 4(d) The largest entanglement is reached at Ji = hi but as J\ diverges from hi, 
C(i,i + 1) decreases. The oscillation appearing at Ji — indicates that the value of hi 
changes the phase of the oscillation. 

There has been great interest in investigating the effect of temperature on the quantum 
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FIG. 5: (Color online) The asymptotic behavior of C(i, i + 1) as a function of (a) A and kT, in units of J\, 
with 7 = 1, ho = hi = 1, and Jq = J\) (b) Ai and kT with 7 = 1, ho = h% = 1, and Jq = 1. 



ent ang lement and the critical behavior of many body systems and particularly spin systems 
(4], 44]- 49]. Osborne and Nielsen have studied the persistence of quantum effects in the 
thermal state of the transverse Ising model as temperature increases 4j. They found that the 
largest amount of entanglement in the system takes place in the region of the parameter space 
close to the critical point A c = 1. This means that at low temperatures the quantum effects 
are still very relevant to the system as manifested by entanglement. Here we investigate the 
persistence of quantum effects under both temperature and time evolution of the system in 
presence of the time-dependent coupling and magnetic field. 



In Fig. 5(a) we reproduce, using our model, the behavior of C(i,i + 1) as a function of 
A and kT for the static case where A = Ao = Ai and ho — hi — 1. As one can see, the 
entanglement is maximum in the vicinity of the critical point A c = 1 and the temperature 
kT = 0. As the temperature increases or A diverges from the critical value, the entanglement 
decays rapidly as the thermal fluctuations destroy the quantum aspects of the system. The 
asymptotic behavior, as t — > 00, of the nearest-neighbor concurrence C(i, i + as a function 



of Ai and kT, while fixing the parameters Jq = 1 and ho — hi — 1, is depicted in Fig. 5(b) 



Interestingly, the entanglement shows a very similar profile to that it manifested in the 
static case, i.e., the system evolves in time preserving its quantum character in the vicinity 
of the critical point and kT = under the time varying coupling. Studying this behavior 
at different values of Jo an d shows that the threshold temperature, at which C(i,i + 1) 
vanishes, increases as Aq increases. Finally, we explore the evolution of next-to-nearest- 



19 



-J„=0.5, J l= 1,kT=0 
J =1, J 1= 1, kT=0 

-J =0.5, J .1, kT.0.1 
J =1, J =1, kT=0.1 



FIG. 6: Dynamics of C(i, i + 2), where t is in units of J 1 , with 7 = 1 and = h\ = 1 for various values 
of J , Ji at fcT = and kT = 0.1. 

neighbor concurrence C(i, i+2), as shown in Fig. [6) As can be noted, C(i,i+2) « C(i,i+1) 
at the same circumstances and vanishes for nonzero temperature. Longer-range concurrence 
C(i, i + r) for r > 3 was found to vanish even at zero temperature. 

V. PARTIALLY ANISOTROPIC XY MODEL 



We now turn to the partially anisotropic system where 7 = 0.5. In this case the x- 
component of the coupling is triple its y component (i.e., J x = 3J y ) and the Hamiltonian 
takes the form 



H 



3J(t) 



N 



i+1 



J(t) N N 



(58) 



i=i 



i=l 



First, we study the dynamics of nearest-neighbor concurrence for this model. In Fig. 7(a) 



we choose the magnetic field to have a constant value of 1 while the coupling parameter is 



0.5 or 2 or a step function changing from 0.5 to 2 (or 2 to 0.5). In Fig. 7(b), we choose 



the coupling parameter to have a constant value of 1 while the magnetic field is 0.5 or 2 or 
a step function changing from 0.5 to 2 (or 2 to 0.5). As one can see, C(i,i + 1) shows a 
nonergodic behavior, similar to the isotropic case, which also follows from the nonergodic 



behavior of the spin correlation functions and magnetization shown in Figs. 7(c) and 7(d) 



Nevertheless, the equilibrium time in this case is much longer than the isotropic case as can 
be seen. 

In Fig. [H]we study C(i,i + 1) as a function of A for different values of J and h and at 
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FIG. 7: Dynamics of the nearest-neighbor concurrence C(i, i + 1), where t is in units of Jj -1 , with 7 = 0.5 
and (a) ho = hi = 1 for various values of Jo and J\ at kT = 0; (b) Jq = J\ = 1 for various values of /io and 
/ii at fcT = 0. (c) Dynamics of the magnetization per spin and (d) the spin-spin correlation function in z 
direction for fixed h = ho = hi = 1 for various values of Jo and J\ at kT = with 7 = 0.5. 

different temperatures. We first study the zero temperature case at different constant values 
of J and h. For this particular case C(i,i+ 1) depends only on the ratio of J and h, similar 



to the isotropic case, rather than their individual values as shown in Fig. 8(a) As can be 
noted, C(i,i + 1) starts from zero, reaches a maximum value at A ~ 0.9, drops to a very 
small value at A ~ 1.1 and then increases rapidly, reaching a constant value for larger values 
of A. The two extremes cases can be explained easily as for h » J, i.e., A << 1, the effect 
of the magnetic field is dominating and the spins are aligned into the z direction and, as a 
result, C(i,i + 1) vanishes. On the other hand, when h « J, i.e., A >> 1, the effect of J 
dominates. However, for this partial anisotropic case, increasing J would increase both J x 
and J y , which causes the spins not to be aligned in a particular direction and consequently 
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FIG. 8: The asymptotic behavior of C(i, i + 1) with 7 = 0.5 as a function of A when ho — hi and Jq = J\ 
at (a) kT = with any combination of constant J and h; (b) fcT = 1 with ho — hi — 0.25, 1,4; (c) kT = 1 
with J = Ji = 0.25, 1, 4; (d) kT = 3 with h a = hi = 0.25, 1, 4. 

C(z,z + 1) maintains an equilibrium finite value. Interestingly, the concurrence shows a 
complicated critical behavior in the vicinity of A = 1, where it reaches a maximum value 
first and immediately drops to a minimum (very small) value before raising again to its 
equilibrium value. The raising of the concurrence from zero as J increases, for A < 1, is 
expected as in that case part of the spins which were originally aligned in the z direction 
change directions into the x and y directions. The sudden drop of the concurrence in the 
vicinity of A = 1, where A is slightly larger than 1, suggests that significant fluctuations 
is taken place and the effect of J x is dominating over both J y and h which aligns most 
of the spins into the x direction, leading to a reduced entanglement value. Studying the 



thermal concurrence in Figs. 8(b), 8(c), and 8(d) we note that the asymptotic value of 



C(i,i+ 1) is not affected as the temperature increases. However, the critical behavior of the 
entanglement in the vicinity of A = 1 changes considerably as the temperature is raised and 
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the other parameters are varied. As illustrated in Fig. 8(b) , the maximum entanglement 
value is reduced and the minimum value reaches zero at high magnetic fields at kT = 1, 
but as the magnetic field is reduced the critical behavior disappears and the entanglement 
makes a direct transition from zero to the equilibrium value where the transition becomes 
sharper and takes place at smaller values of A as we increase the magnetic field. A similar 



behavior is shown in Fig. 8(c) where for small values of the coupling J the critical behavior 



disappears as well. The effect of higher temperature is shown in Fig. 8(d) where the critical 



behavior of the entanglement disappears completely at all values of h and J, which confirms 
that the thermal excitations destroy the critical behavior due to suppression of quantum 
effects. 























FIG. 9: (Color online) C(i,i + 1) as a function of Ai and t, in units of J^ 1 , at kT = with 7 = 0.5, 
ho = hi = 1 and (a) Jo = 1; (b) Jo = 5. The asymptotic behavior of C(i, i + 1) as a function of (c) Jo and 
Ji at kT — with 7 = 0.5 and ho = h\ = 1; (d) ho and hi at kT = with 7 = 0.5 and Jo = Ji = 1, where 
ho, h%, and Jq are in units of Jj. 
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In Fig. [9] we investigate the dependence of the time evolution and asymptotic behavior 
of concurrence on the different magnetic field and coupling parameters. We have studied 
C(i,i + 1) as a function of Ai and t for many different selections of h = ho = h% and Jq and 
found that the concurrence behavior depends mainly on whether Aq > 1 or < 1. To test 



that behavior, Figs. 9(a) and 9(b) show C(i,i + 1) as a function of Ai and t with J = 1 and 



5 respectively for fixed ho = hi = 1 and kT = 0. As can be seen in Fig. 9(a), where Jo = 1 
(Ao = 1), the concurrence at any time t > starts with a finite value (or zero) at Ai = Ao and 
then decays to zero but increases again, reaching a maximum value in the vicinity of Ai = 1 
and finally vanishes permanently as Ai increases. Interestingly at Ai = the concurrence 
shows an oscillatory behavior in time, as was the case in the completely anisotropic model. 
For values of Ai around the critical value, the concurrence approximately maintains its initial 
value as time elapses. However, for larger values of Ai, the concurrence starts initially with 
a finite value but decays sharply to zero in a very short period of time. On the other hand, 



in Fig. 9(b) where we set Jo = 5 (i.e., Ao = 5), at any time t > the concurrence starts 
with a finite value at Ai = and increases rapidly as Ai increases because increasing J\ 
reduces the alignment in the z direction. The concurrence C(i,i + 1) reaches a maximum 
value at Ai ~ 0.2 and vanishes at Ai ~ 1. Finally, the concurrence increases, reaching a 
constant value for Ai ~ 2 or larger. The variation of the concurrence in this case with time 
is very limited, as one can see; it approximately maintains its initial value especially for all 
Ai < 1 and Ai > 3. This critical dependence of the concurrence dynamics on the initial 



value of the coupling parameter is emphasized in Fig. 9(c), where the asymptotic value of 



the concurrence is depicted as a function of both J and J\. In Fig. 9(d) the asymptotic 
behavior of the concurrence is explored as a function of ho and hi while fixing the coupling 
as J\ — Jo — 1. The behavior of the concurrence is very close to the completely anisotropic 
case except for the region where < ho < 1 and < h\ < 1 where the concurrence starts 
with a finite value at ho = hi = and decays gradually until it vanishes at ho = hi « 1. 
For higher values of the coupling J\ = Jq the rate of decay of the concurrence every where 
is smaller and the peaks are broadened. 

The persistence of quantum effects as temperature increases and time elapses in the 



partially anisotropic case is examined and presented in Fig.[10j In Fig. 10(a) the concurrence 
is plotted as a function of A = Ao = Ai and kT with ho — hi — 1 and Jo = J\. As one 
can see, the concurrence shows the expected behavior as a function of A and decays as the 



24 



FIG. 10: (Color online) The asymptotic behavior of C(i,i + 1) as a function of (a) A and kT, in units of 
Ji, with 7 = 0.5, ho = hi = 1, and Jo = Ji; (b) Ai and fcT with 7 = 0.5, /to = fti = 1, and Jo = 1. 



temperature increases. As one can see, the threshold temperature where the concurrence 



vanishes is determined by the value of A, it increases as A increases. In Fig. 10(b) the 
asymptotic behavior of the concurrence as a function Ai and kT is illustrated. Clearly the 
nonzero concurrence shows up at small values of kT and X\. The concurrence has two peaks 
versus Ai but as the temperature increases, the second peak disappears. Very interestingly, 
the first peak raises as temperature increases then decays again and vanishes for kT m 0.9 
or larger. 
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FIG. 11: Dynamics of C(i,i + 2), where t is in units of J 1 , with 7 = 0.5 and ho = hi = 1 for various 
values of J , Ji at kT = and kT = 0.25. 

The behavior of the next-to-nearest-neighbor concurrence C(i, i + 2), is shown in Fig. [11] 
As expected, C(i, i + 2) << C(i,i + 1) at the same circumstances. Studying longer-range 
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concurrence C(i,i + r) shows that it vanishes for r > 3. 

VI. ISOTROPIC XY MODEL 

In this section we consider the isotropic system with 7 = (i.e., J x = J y ); in this case 
the Hamiltonian assumes the form 

JV 



H 



4r E °f - E ^ - E WW 

/ i=l / i=l i=l 



(59) 



We start with the dynamics of the nearest-neighbor concurrence; in Fig. [T21 we first 
choose the magnetic field to have a constant value of 1 while the coupling parameter is 2, or 
a step function changing between 0.5 and 2. We also study the case with a constant coupling 
parameter of 1 while the magnetic field is 2, or a step function changing between 0.5 and 
2. Testing the concurrence for several different values of the magnetic field and coupling 
parameter, we note that C(i, i + 1) takes a constant value that does not depend on the final 
value of the coupling J\ and magnetic field h\. This follows from the dependence of the spin 



correlation functions and the magnetization on the initial state only as shown in Fig. 12(b) 



This is because the initial coupling parameters J x and J y , which are equal, force the spins to 
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FIG. 12: Dynamics of the nearest- neighbor concurrence C(i,i + 1), where t is in units of J[~ , with 7 = 
for various values of Jo,Ji,ho, and hi at kT = 0. (b) Dynamics of the magnetization per spin and the 
spin-spin correlation function in the z direction for fixed h = ho = hi = 1 for various values of Jo and Ji at 
kT = with 7 = 0. 

be equally aligned into the x and y directions, apart from those in the ^-direction, causing 
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a finite concurrence. Increasing the coupling parameters strength would not change that 
distribution or the associated concurrence at constant magnetic field. 
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FIG. 13: The asymptotic behavior of C(i, i + 1) with 7 = as a function of A when ho = hi and Jo = Ji 
at (a) kT = with any combination of constant J and h; (b) kT = 1 with ho = hi = 0.25, 1, 4; (c) kT = 1 
with Jo = Ji = 0.25, 1,4; (d) kT = 3 with h = hi = 0.25, 1,4. 



In Fig. [13] we study C(i,i + 1) as a function of time-independent A for different values 
of J = Jq = J\ and h = ho = h\ and at different temperatures. Again C(i,i + 1) depends 
only on the ratio of J and h rather than their individual values at kT = as shown in 



Fig. 13(a) As can be seen, C(i, i + 1) starts from zero, increases rapidly at A = 1, and then 
maintains a constant value as A increases. In this case, when h » J, i.e., A << 1, the 
effect of magnetic field dominates, causing the spins to be aligned to the z direction and 
as a result C(i,i + 1) vanishes. On the other hand, when h << J, i.e., A >> 1, the effect 
of the coupling dominates and the spins are equally aligned in both x and y directions and 
the concurrence maintains a constant finite value. Interestingly, raising the temperature as 
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shown in Fig. 13(b) does not reduce the concurrence but causes the values of the coupling 
and the magnetic field to affect the concurrence independently at nonzero temperature. Also 



in Fig. 13(b) , one can notice that decreasing h at kT = 1 causes the change in C(i,i + 1) to 



be less rapid and to reach equilibrium value at larger values of A. In Fig. 13(c), we fix the 



coupling at different values and study the concurrence as the magnetic field changes. The 
concurrence C(i,i + 1) vanishes when J < 1, i.e., when the magnetic field dominates, while 
for J > 1 it manifests the same behavior as before. The effect of higher temperatures is 



shown in Fig. 13(d) where it causes smoother change in the concurrence. 



(a) 









FIG. 14: (Color online) (a) C(i,i + 1) as a function of Ai and t, in units of J 1 , at kT — with 7 = 0, 
ho = hi = 1 and Jo = 5; (b) C(i, i + 1) as a function of Ao and t at kT = with 7 = 0, ho = hi = 1 and 
Ji = 5. 

The time evolution of nearest-neighbor concurrence as a function of the time-dependent 



coupling is explored in Fig. [T3J where we fix the magnetic field. Figure 14(a) shows C(i, z + 1) 
as a function of Ai and t where ho — h± — 1 and Jo = 5 at kT = 0. Clearly, C(i,i + 1) 
is independent of Ai. Studying C(i,i + 1) as a function of A and t with h — hi = 1 at 



kT = for various values of Ji, we note that the results are independent of J\. Figure 14(b) 



represents the case where Ji = 5. Again, as can be noticed when Jo < ho, the magnetic 
field dominates and C(i, i + 1) vanishes. While for J > h , C(i,i+ 1) has a finite value, as 
discussed above. 



The effect of temperature on concurrence is investigated in Fig. [151 In Fig. 15(a) we plot 
the asymptotic concurrence C(i,i+ 1) as a function of Ao and kT. Clearly, as kT increases 
the threshold Ao, at which C(i, i + 1) starts to have a finite value, increases. An interesting 
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FIG. 15: (Color online) The asymptotic behavior of C(i, i + 1) as a function of (a) Ao and kT, in units of 
Ji, with 7 = 0, ho = hi = 1, and J\ — 1; (b) kT with 7 = 0, ho = hi = 1, and Jo = Ji = 1. 

behavior of the concurrence is featured here, studying the asymptotic concurrence C(z, i + 
as a function of Ai and kT, one observes that for A < 1, C(i, i + 1) starts from at /cT = 0, 
grows up as kT increases, reaching a maximum value at kT 0.3, and then vanishes again 



for kT w 0.9 as shown in Fig. 15(b) In Fig. [161 we depict the time evolution of next-to- 
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FIG. 16: Dynamics of C(i, i + 2), where t is in units of J Y , with 7 = and /lo = ^i = 1 f° r various values 
of Jo, Ji at fcJ = and kT = 0.1. 



nearest-neighbor concurrence C(i, i + 2). As expected, C(i,i + 2) « C(i,i + 1) at the same 
circumstances. The longer-range concurrence C(i,i + r) vanishes for r > 3. 

Finally, we explore the asymptotic behavior of the nearest-neighbor and next-to-nearest- 
neighbor concurrence in the A-7 phase space of the one-dimensional XY spin system under 



the effect of a time-dependent coupling J(t). In Figs. 17(a) and 17(b) we plot C(i,i + 1) 
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and C(i,i + 2) respectively as a function of the parameter Ai and the degree of anisotropy 
7 for constant magnetic field ho — hi — 1 and Jo — 1 at kT = 0. As one can notice, the 
nonvanishing concurrences appear in the vicinity of A = 1 or lower and vanishes for higher 
values. One interesting feature is that the maximum achievable nearest-neighbor concurrence 
takes place at 7 = 1, i.e., in a completely anisotropic system, while the maximum next-to- 
nearest-neighbor concurrence is achievable in a partially anisotropic system, where 7 ~ 0.3. 




FIG. 17: (a) (Color online) The asymptotic behavior of C(i, i + as a function of Ai and 7 with ho = h\ = 1 
and Jo = 1 at kT = 0. (b) The asymptotic behavior of C(i,i + 2) as a function of Ai and 7 with ho = h\ = 1 
and J = 1 at kT = 0. 

VII. CONCLUSIONS AND FUTURE DIRECTIONS 

We have investigated the entanglement evolution in an infinite one-dimensional XY model 
in an external time-dependent magnetic field at zero and finite temperature. The nearest- 
neighbor interaction between the spins were considered time dependent. An exact solution 
was presented for a step function form of both the time- dependent coupling and magnetic 
field. The system showed nonergodic and critical behavior at all degrees of anisotropy and 
for all different choices of coupling and magnetic field. At zero temperature and constant 
magnetic field and coupling, the asymptotic behavior of the system at the infinite time limit 
depends only on the ratio of the coupling to the magnetic field, not their individual values 
but changes for nonzero temperature. For many system setups the initial values of the 
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coupling and magnetic field dictate the asymptotic behavior of the entanglement regardless 
of their final values at different degrees of anisotropy and may lead to asymptotic residual 
entanglement in the system. Interestingly, studying the dynamics of entanglement at zero 
and finite temperature showed that the quantum properties of the system are preserved 
within certain regions of the coupling, magnetic field, and temperature space that vary 
significantly depending on the degree of anisotropy of the system. Particularly, the quantum 
effects in the transverse Ising model persist in the vicinity of both its critical phase transition 
point and zero temperature as it evolves in time. In future work, it would be interesting to 
study the XY model under the effect of time- dependent coupling and magnetic field where 
the function form for each one could differ from the other and take other forms of practical 
interest, such as the sinusoidal and exponential. This would need in that case the use of 
numerical methods along with the analytical ones to treat the system. 
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